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Abstract 

A new approach is proposed, namely CSSF MIMO radar, which applies the technique of step 
frequency (SF) to compressive sensing (CS) based multi-input multi-output (MIMO) radar. The proposed 
approach enables high resolution range, angle and Doppler estimation, while transmitting narrowband 
pulses. The problem of joint angle-Doppler-range estimation is first formulated to fit the CS framework, 
i.e., as an t\ optimization problem. Direct solution of this problem entails high complexity as it employs 
a basis matrix whose construction requires discretization of the angle-Doppler-range space. Since high 
resolution requires fine space discretization, the complexity of joint range, angle and Doppler estimation 
can be prohibitively high. For the case of slowly moving targets, a technique is proposed that achieves 
significant complexity reduction by successively estimating angle-range and Doppler in a decoupled 
fashion and by employing initial estimates obtained via matched filtering to further reduce the space that 
needs to be digitized. Numerical results show that the combination of CS and SF results in a MIMO 
radar system that has superior resolution and requires far less data as compared to a system that uses a 
matched filter with SF. 

Keywords: Compressive sensing, MIMO radar, step frequency, DOA estimation, Doppler estimation, 
range estimation 

I. Introduction 

Benefiting from the idea of multiple-input multiple-output (MIMO) communication systems, MIMO 
radar systems have received considerable attention in recent years. A MIMO radar [l]-[3] although having 
some similarities with phased-array radar, is different in that it transmits multiple waveforms from its 
antennas, and in general offers a more flexible configuration. A MIMO radar with widely separated 
antennas [4] views a target from different decorrelated directions and thus enjoys spatial diversity. For 

1 This work was supported by the Office of Naval Research under Grant ONR-N-00014-07-1-0500, and the National Science 
Foundation under Grants CNS-09-05398 and CNS-04-35052. 
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MIMO radar with colocated antennas [5] [6], independent waveforms enable superior angular resolution 
as compared to phased-array radar. We focus on the colocated type of MIMO radar throughout this paper. 

Compressive sensing (CS) is a relatively recent development [7]-[10] and has already been applied 
successfully in diverse fields such as image processing and wireless communications. The theory of CS 
states that a A'-sparse signal vector x of length N can be recovered exactly with high probability based 
on 0(K log N) measurements via i\ -optimization. Let *P denote the basis matrix that spans the space in 
which the signal is sparse, and O be the measurement matrix that is used to linearly compress the signal. 
The 1 1 -optimization problem is formulated as follows: 

min||s||i, s.t. to y = Ox = OTs (1) 

where s is an Nx 1 sparse vector with K principal elements while the remaining elements can be ignored; 
$ is an M x W matrix with M <^ N. The product 0*P is referred to as the sensing matrix. The uniform 
uncertainty principle (UUP) [9] [11] indicates that if every set of sensing matrix columns with cardinality 
less than the sparsity of the signal of interest is approximately orthogonal, then the sparse signal can be 
exactly recovered with high probability. This implies that O is incoherent with *P. Throughput this paper, 
we will refer to the solution of a problem of the form (1) as the CS approach. 

By exploiting the sparsity of radar signals in various spaces, CS has been applied to radar systems 
[12]-[15] and to MIMO radar [16]-[19]. In [17] and [18], a uniform linear array was considered as a 
transmit and receive antenna configuration and the CS approach was applied using a submatrix of the 
identity mattix as the measurement matrix. Bounds on the achievable range and azimuth resolution and 
the number of recoverable targets were derived in [18]. In [16] and [19], the authors proposed a CS-based 
MIMO radar system implemented on a small scale network. According to [19], randomly located network 
nodes, each equipped with a single antenna, serve as transmit and receive antenna elements. The transmit 
nodes transmit periodic pulses. The receive nodes forward their compressively obtained measurements to 
a fusion center. Exploiting the sparsity of targets in the angle-Doppler space, an i\ -optimization problem 
is formulated and solved at the fusion center in order to extract target angle and Doppler information. 
This approach achieves the superior resolution of MIMO radar with far fewer samples than required by 
conventional approaches, which implies lower power consumption for the receive nodes. 

The application of CS to step-frequency radar (SFR) [24] was investigated in [25]-[27]. SFR transmits 
pulse trains of varying frequency. Thus, although the individual pulses are relatively long in duration 
and are narrowband, the transmit signal is effectively wideband. Since range resolution increases with 
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the signal bandwidth, SFR achieves high range resolution. At the same time, SFR does not suffer from 
low signal-to-noise ratio (SNR) at the receiver typically associated with wideband systems that rely on 
short duration pulses. In [27], it was found that the CS approach can significantly reduce the number of 
pulses required by SFR to achieve a certain resolution. A CS -based data acquisition and imaging method 
was proposed in [25] for stepped-frequency continuous-wave ground penetrating radars, and in [26] CS 
was applied to stepped-frequency through-the-wall radar imaging. In both cases it was shown that the 
CS approach can provide a high-quality radar image using many fewer data samples than conventional 
methods. 

In this paper, we consider a more general scenario than that of [16] and [19]. The methods of [16] 
and [19] assume that the targets are located in a small range bin and the sampling is synchronized with 
the first target return. Such assumptions do not allow for range estimation. In this paper, the targets 
can be located across several range bins. We propose CSSF MIMO radar, an approach that applies step 
frequency to CS-based MIMO radar. Two types of CSSF MIMO radar systems are considered, i.e., 
linear step-frequency radar (LSFR), and random step-frequency radar (RSFR), and their effects on the 
CS approach are studied. The proposed approach enables high resolution range as well as angle and 
Doppler estimation. The problem of joint angle-Doppler-range estimation is first formulated to fit the 
CS framework, i.e., as an €\ optimization problem. Solving this problem entails high complexity as 
it employs a basis matrix whose construction requires discretization of the angle-Doppler-range space. 
The complexity increases with the size of the basis matrix, or equivalently, as the discretization step 
decreases; the latter step needs to be as small as possible as it sets the lower limit of resolution. For 
slowly moving targets, a technique is proposed that successively estimates angle -range and Doppler in a 
decoupled fashion, and employs initial estimates obtained via a matched filter (MF) to further reduce the 
space that needs to be digitized. In [20] and [21], information on the support of a sparse signal was used 
in the minimization process resulting in complexity reduction. In our case, we do not explore the role of 
initial estimates on the minimization process, as they are not expected to be very reliable. Instead, we use 
them only as a guide for the construction of the basis matrix. The preliminary results of CSSF MIMO 
radar and the decoupled scheme were published in [22] and [23] which consider the same signal models 
as in [16] and [19]. This paper extends the work of [22] and [23] to the general scenario aforementioned 
and offers an mathematic insight into CSSF MIMO. 

The paper is organized as follows. In Section II we provide the signal model of a CS-based MIMO 
radar system. In Section III, we introduce the proposed CSSF MIMO radar system. A decoupled scheme 
for CSSF MIMO is described in Section IV. Simulation results are given in Section V for the case of 
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slowly moving targets. Finally, we make some concluding remarks in Section VI. 

Notation: Lower case and capital letters in bold denote respectively vectors and matrices. The expecta- 
tion of a random variable is denoted by £{•}. Superscripts {-) H and Tr(-) denote respectively the Hermitian 
transpose and trace of a matrix. A(m, n) represents the (m, n)th entry of the matrix A. 

II. Signal Model for CS-based MIMO Radar 

Let us consider the same setting as in [19]. In particular, assume K point targets and colocated transmit 
(TX)/receive (RX) antennas that are randomly distributed over a small area. The k-th target is at azimuth 
angle 6 k and moves with constant radial speed v k - Let M, a'-)/^, o-p denote the location of the i-th TX/RX 
node in polar coordinates. The number of TX and RX nodes is denoted by M, and N r , respectively. Let 
d k {t) denote the range of the k-th target at time t. Under the far-field assumption, i.e., d k {t) » r'/ r , the 
distance between the /th transmit/receive node and the k-th target d' ik /d'- k can be approximated as 

d'/ k r (t) a d k (t) - rfl r m = 4(0) - rf!'\e k ) - v k t (2) 

where rfj r {0 k ) = r/ r cos(% - af). 

Assuming that there is no clutter, each TX node transmits periodic pulses of duration T p and pulse 
repetition interval (PRI) T. The target return from the k-th target arriving at the /-th antenna during the 
m-th pulse is 

M, 

y]jfi = J] fikX.it - (4(0 + d' lk {t))lc) exp(j2nf(t - (d' lk (t) + d r lk (t))/c)) (3) 
(■=] 

where c, f and fi k denotes the speed of light, the carrier frequency, and the reflection coefficient of 
the fc-th target, respectively; Xj(t) represents the transmit waveform of the i-th node. The demodulated 
baseband signal corresponding to a single target can be approximated by 

M, 

vL(f) x Z PkXi{t ~ 2fi W°)/ c ) exp(- 7 -27r/(4(0 + d r lk (t))/c). (4) 
(■=] 

In the above equation, the time delays in the received waveforms due to the £-fh target are all the 
same and equal to 2d k (0)/c. This approximation is enabled by the assumption of narrowband transmit 
waveforms, slowly moving targets and colocated nodes. The fact that the targets can be in different range 
bins implies that the delays corresponding to different targets will be different. Therefore, sampling the 
received signal and ensuring that all target returns fall in the sampling window would require time delay 
estimation. However, in a CS scenario, there are not enough data to obtain such estimates using traditional 
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methods, e.g., the matched filtering method (MFM). In the following, we will extract the range and angle 
information simultaneously using the CS approach without assuming availability of time delay estimates. 
The compressed samples collected by the Z-th antenna during the m-th pulse can be expressed as 

K 

rim = Yj he i2np '"' k ^il)(Jk)C n Xy{9 k ) + <t> t n lm (5) 

where 

1) Pimk = ~ 2d f )f + n -^~+fk{m - l)T, where f k = ^fis the Doppler shift induced by the jt-th target; 

diag{X ff X} = [1 l] r ; IT S ,1 = 0, . . . ,L - 1, represent the time within the pulse (fast time) and 

thus the pulse duration is T p = LT S ; 

2) O/ is the M x (L + L) measurement matrix for the 1-th receive node where LT S is the maximum 
time delay and known in advance. The measurement matrix has elements that are independent and 
identically distributed (i.i.d) Gaussian random variables; 

3) \(8 k ) = [ e j^^\ ...,e'^< m f and D(f k ) = diag{[e^ or \ . . . ,e j2Kfk(L - 1)T ']}; 

4) r k = L^f^J and C Ti = [0z. xri , I/., Lx( £_ n) ] r . Here, we assume that the target returns completely fall 
within the sampling window of length (L + L)T S , and that T s is small enough so that the rounding 
error in the delay is small, i.e., xi(t - r k ) w xi(t - L 2 ^ (Q) J). 

5) n/ m is the interference at the l-th receiver during the m-th pulse, which includes a jammer's signal 
and thermal noise. 

Let us discretize the angle, speed and range space on a fine grid, i.e., respectively, [Si, . . . 
[bi, . . . ,b Nb ] and [c\, . . . , cjvJ. Let the grid points be arranged first angle-wise, then range-wise, and 
finally speed-wise to yield the grid points (a n ,b n ,c n ),n = l,...,N a NbN c . Through this ordering, the grid 
point (a na ,b nb ,c nc ) is mapped to point (a n , b n , c n ) with n = (n^ - \)n a n c + (n c - \)n a + n a . We assume that 
the discretization step is small enough so that each target falls on some angle-speed-range grid point. 
Then (5) can be rewritten as 

rim - <S>i [j] s n eflm " D f^T") C L ^jXv(a„) + 

\n=l ' ' 



(6) 



reflection coefficient of the target, if there is a target at (a n ,b n ,c n ) 
where s n = \ ,N = N a N h N c , and 

0, if there is no target at (a n ,b n ,c n ) 

-2c„f rfMn)f 2b n f(m - \)T 

qimn ~ + + ■ (7) 
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In matrix form we have 

r hn = Zm s + 0,n /m (8) 

where s = [s\, s^] T and 

©/m - ®i [e^Yy{2b,flc)C v ^y{a{), e ^"» w D(2^// C )C L ^ J Xv(a JV )] . (9) 

According to the CS formulation, 0/ m is the sensing matrix and *P/ m is the basis matrix. 

If the number of targets is small as compared to N, the positions of the targets are sparse in the 
angle-speed-range space, i.e., s is a sparse vector. The locations of the non-zero elements of s provide 
information on target angle, speed and range. 

All the receive nodes forward their compressed measurements to a fusion center. We assume that the 
fusion center has the ability to separate the data of different nodes from each other. This can be done, for 
instance, if the nodes send their data over different earners. The fusion center combines the compressively 
sampled signals due to N p pulses obtained at N r receive nodes to form the vector r. It holds that 

r = [r T n ,...,r T 1Np ,...,r T NrN f = 0s + n (10) 

where = [(0 n ) r , . . . , (& lNp ) T , (Q Nl . Np ) T ] T and n = [(<D inil ) r , . . . , (Qm 1JVp ) r , . . . , n NrNp ) T ] T . 

Subsequently, using the predefined measurement matrices, O/, 1=1, ...,Af,-, based on the discretization 
of the angle-speed-range space, and also based on knowledge of the waveform matrix X, the fusion center 
recovers s by applying the Dantzig selector [28] to (10) as 

§ = min||s||i s.t. ||0 H (r - 0s)|U < V- (11) 

According to [28], the sparse vector s can be recovered with very high probability if fi = (l+t~ l ) -\/21og N 'a 1 cr, 
where t is a positive scalar, cr max is the maximum norm of columns in the sensing matrix 0, and cr 2 is 
the variance of the interference in (10). A numerical method to determine the value of /i is described in 
[28]. 

III. Introducing Step Frequency to CS-MIMO radar 
Let us consider a MIMO radar system in which the carrier frequency of the m-th pulse equals 

f,n=f + Hn (12) 
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where / is the center carrier frequency and Af m denotes the frequency step, m - 1, . . . , N p . 
The baseband samples collected by the l-th antenna during the m-th pulse can be expressed as 

K 



rim = <S>iY J pke j2npM B(fmk)C Tk Xv m (9 k ) + <D,n /m (13) 



k=\ 

where 

c 

-2d k {0)f m , rf,mf m , . . nT n .. 

and pi^ = + + f mk (m - V)T. (14) 

c c 

Then, based on discrete grid points of the angle-speed-range space, (13) can be rewritten as 

f lm = <D,T /m s + 0,n /m 

= 0, m s + 0/n /m (15) 

where 

T /ffl - [e' 2 ' r *"'D(2^ 1 / m /c)C L ^ J Xv m (a 1 ), . . . , ^'»'"D(2fe iV / m /c)C L ^ J Xv m (fl Af )], 
-2c n f m rf^Qfn 2b n f m (m - l)T 

C C C 

and & lm - O/T,,,,. (16) 

At the fusion center, the compressively sampled signals due to N p pulses obtained at N r receive nodes 
are stacked as 

f =0s + n (17) 

where 

= K@ n ) T , (&wf, (& Nr N,f] T . (18) 

Recovery of s is performed as in (11) where is replaced with 0. 
In the remainder of the paper, we make the two assumptions: 

• (Al) The targets are slowly moving. Therefore, the Doppler shift within a pulse can be ignored, i.e., 

f m (2T s (L + L- l)b n )/c * 0, n = 1, . . . ,N. 
. (A2) The radar waveforms are independent across transmit nodes and thus J Xj(t)x*,(t + T)dt, i + i 



January 17, 2011 



DRAFT 



is negligible as compared to J _ Q + r)Jf. 

A. Range resolution 

In this subsection we study the relationship between range resolution and the ambiguity function. For 
the conventional radar systems that uses a matched filter to extract target information, the ambiguity 
function (AF) characterizes the response to a point target and determines resolution. Let us assume that 
there is a target at (9,d,v). The matched filter looking for a target at {ff ,d' ,v') yields 

^JL ML E]L ^ w) ^r m _,f Ae , ) _ l] r {f) , ) _ 2M 

X (Ad, Av, 6, 0) = J] Yj JjXisAM, Av)e }2nf " e (19) 

1=1 ,/=l ,„=1 

where Ad = d - d', Av = v - V and 

Xu'AAd, Av) = Jxi(t)x*,(t + 2Adlc)e j2nf '" 2 -^'dt. (20) 

Equation (19) is the AF for SF MIMO radar, where SF MIMO radar refers to MIMO radar that uses the 
SF technique. Unlike the AF for MIMO radar [29], the carrier frequency is varying between pulses in 
(19). 

To investigate the range resolution let us set Av = and 9 = 0'. Then, the AF becomes 

M, N„ 

x(Ad,o,e,e) = N r J] 2_jX^ M 'W e — 1 — 

i,i'=l m=l 

= N r 2^ e iW«(-2Ad/c) ^ Xi (f)x*,(t + 2 Ad 10 dt + N r 2_j 2_j e j2nfn '~ Xi(t)x*,(t + 2Ad/c) dt 

m=\ i=i' J' m=\ i+i' ^' 

(21) 

Due to (A2), the term Ax{Ad) is negligible as compared to the product x\(Ad)x2{ Ad) in (21). One can 
see that^i(AJ) and xi(Ad) are respectively the AF of SF single-input single-output (SISO) radar and 
MIMO radar, both for Av = and 6 = 6'. It can seen from (21) that a colocated MIMO radar has no 
gain on range resolution as compared to a SISO radar, i.e., the range resolution of MFSF MIMO radar 
is at least equal to the best between the range resolution of SF SISO radar and SISO radar, where MFSF 
MIMO radar refers to matched filter based MIMO radar that uses the SF technique. 

In [30], in a study of CS -based SISO radar, it was observed that the maximum value of the correlation 
of two different columns of the basis matrix is equal to the second largest value of the discrete AF surface. 
The recovery performance of CS approaches, however, is directly related to the column correlation of 
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the sensing matrix rather than the basis matrix. Unlike [30], we next study the relation of the AF and 
the column correlation of the sensing matrix for the proposed CSSF MIMO radar. This analysis will 
provide a clue for comparing the resolution of CS and matched filter in the context of SF MIMO radar, 
i.e., CSSF MIMO radar and MFSF MIMO radar. 

On letting p^, denote the column of the sensing matrix corresponding to the k-th grid point in the 
angle-speed-range space, we have 

< p*,p*< > = £fy 2 ^" 

/=1 m—\ cTs * ' ~£ ' ' ^ ' / 

g* 

N, N p L+L 

= ZZ Z ^^^sMgMMP,*) 

/=1 m=l p,q=l 

N r N P L+L M, 

1=1 m=l p,q=l i,i'=l 

■ x, L - l)T, - —) 4 \(P ~ - — ) (22) 



c l \ c 

where Ad kk > = c k _ ct and Av^ = b k - b k <. For simplicity, in the above we assumed that the receive 
nodes use the same measurement matrix; thus the index / was dropped in A. 

Taking the elements of the measurement matrix <D to be independent and Gaussian N(0, jj-j), the 
expectation of < p^, p k > > with respect to the elements of <D equals 

m N ' N " M < 
E{<p k ,p k > >} = Z Z Z e j2nf "' (n '- (ak)+ ^ (ak) ~'''i' (a ^ 

L+L 

• s ( P -mv W )/c 



ocx(Ad kk >,Av kk >,a k ,a k >). (23) 

One can see from the above equation that the expectation of the column correlation of the sensing matrix is 
proportional to the discrete AF. To focus on the range resolution we set a k = a k > and Av kk , = 0. Essentially, 
the range resolution of MFSF MIMO radar corresponds to the smallest range difference between two 
targets, Ad kk >, that sets the AF to zero. Based on the UUP in [11], however, the coherence of the sensing 
matrix does not have to be zero for exact recovery; a small level of coherence is good enough. Therefore, 
CS -based radar systems have the potential to improve range resolution. This possibility will be confirmed 
via simulations in Section V (see Fig. 2). 
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B. The effect of signal bandwidth on CSSF-MIMO radar 

In an LSFR system, the carrier frequency increases by a constant step between pulses, i.e., Af m = 
(m - l)Af. This type of SF radar can be efficiently implemented using the Inverse Discrete Fourier 
Transform (IDFT) [24]; however, it suffers from range ambiguity if the distance between a target and 
receive nodes exceeds the value R u - The range ambiguity can be removed by randomly choosing 
the step frequency within a fixed bandwidth at the expense of increased sidelobe as compared to the 
LSFR [32]. In this section, we investigate the effect of the number of pulses N p (or equivalently, the 
bandwidth) on range resolution for two types of CSSF MIMO radar, i.e., LSFR and RSFR, in terms of 
the coherence of the sensing matrix (see (29)). Consistent with [32], which discussed convectional radar 
systems using the MFM, we find that the RSFR requires more pulses than LSFR to achieve the same 
range resolution for CS -based MIMO radar. 

Since an increase in the number of receive nodes does not improve the range resolution, for simplicity 
we consider one receive node only. The correlation of columns p/ f and for a* = ay and = by equals 



Pkt 



N p L+L 

ZZ« 



m=l p,q— 1 
L+L 



( M, M, 

Qkk> (m, P, q, i, i) + 2j Qkk' (m, p, q, i, i') 



^ e j2nU-2Ad^)/c ^ A(p,q)e j2r!fl " i2Tshk(q ~ p)yc 



in— 1 



/=1 



(24) 



where 



Qkk<m, p, q, i, /') = e V«m^-rtM^ Xi [i^q - l)T s - ** L - l)T s - ^ J . (25) 
Due to (Al) and the discretized version of (A2), we can ignore the Doppler shift within a pulse and 



the second term Qkk'(m, p,q,i,i') in (24). Therefore, (24) becomes 



Pkt 



N p L+L M, , v , 

£ e^^c £ A(p, q) J] Xi L -Ws-^U Up ~ »T, - ^ 



(26) 



Eq. (26) can be rewritten as 



Pkt 



NpPkk 



m- 1 



j a kk' (f+Afm) 



MM) 



k = k' 
Pkk> k±k' 



(27) 
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where Af = [Af u . . . , Af Np ], 



Pkk> 



L+L M, 



and a kk > = _ 



£ Y j A{p,q)x i {{q - l)T s - -^)x*((p - l)T s - -^) 

p,q=i i- 1 

AnAd kk , 



(28) 



Then, the coherence of the sensing matrix corresponding to columns and can be written as 

p kk > h{Ai)p kk > 



Pkk>(®) = 



(29) 



y/PkkPk'k' N p y/pkkPk'k' ' 

1) Linear step frequency: If the carrier frequency increases by a constant step Af between adjacent 
pulses, i.e., Af m = (m - 1)A/, then 



Pkk>(@) 



|1 _ e ia w ^v\p kk , | sinC^A/Ayi 



N„ 



(30) 



|1 _ e j ai *' A f\N p -^pidpTiP »p 
It can be easily seen that an increase in N p tends to reduce the coherence and thus improves the range 
resolution. 

Let arf?. ., denote the travel- time difference between the signals sent from the transmit node i to the 

kk'u' ° 

target located at the kth grid point at time instant pT s , and from the transmit node i' to the target located 
at the k'th grid point at time instant qT s . It holds that 



a p k l w = {-2Ad w + 2T s b k (q ~ p) + q\(a k ) - q'M^/c. 



Regarding the approximation error, the term discarded in (26) is 



(31) 



L+L 



M, 



Pkk< = Y j e^~ 2Ad ^ 2 ^q^^^^Que^p^Uf) 

p,q=l i+i' 



m-\ 
L+L 



= L " kk " Mp, q) ^ Xi (* - 1)7, - - x, (p - l)r. - — . (32) 

The amplitude of e . 2jA " e } T is bounded by ATp. For independent waveforms, the approximation 
error p kk! in (24) is always negligible as compared to p kk >. 

Let p t denote the maximum coherence of that guarantees exact recovery of the sparse vector with 
high probability via the Dantzig selector. The minimum number of pulses required to achieve a certain 
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resolution can be obtained by solving 



N* p = min N p 



s.t. 



|1 - e** A/ "'|pj* 



\l-e*>**f\N p jpZpZi 
k,k' = l,...,N aadk±V. 



<P-t 



(33) 



The above problem is easy to solve, for example by trying different values for N p ; however, it requires 
a value for p. t . In [18], a rough estimate of //, in the presence of mild interference was offered. In general, 
H, must be determined experimentally 

2) Random step frequency: Assuming that the frequency steps over pulses are i.i.d uniform random 
variables, i.e., A/„, ~ U(0,2b), the expectation of square coherence over Af m is given by 

2 



£{/4(0)} = £ 



2^1 e*"<f+*M Pa , 



Put 



PkkPk'k' 



NjpkkPk'k' 

1 N 2 p -N pS in 2 (a w b)) 



N„ 



N 2 
p 



< b2 



(34) 



For a fair comparison, we set LSFR and RSFR to cover the same frequency band, i.e., set b equal to 
Af(N p - l)/2. Then (34) can be rewritten as 



£{/4(©)} 



PkkPk'k' N p 

P\k' 
PkkPk'k' N p 



1 + 



1 + 



' 2 



kN p -l)akk'&f) \ 



sin 2 (i(Af p -l)tr H . 



{N p -l){2nAfhd kk 



'A/) \ 



(35) 



As the term (2nAfAdtk'/c) 2 increases, the expected value of the squared coherence becomes approxi- 
mately equal to l/N p . This holds when the product of radian frequency step 2nAf and the range spacing 
of grid points A<4^ is comparable to the speed of light c. 

Since the coherence of the sensing matrix for RSFR cannot be obtained directly, we instead compare 
the squared coherence of the sensing matrix for LSFR and RSFR. For large N p , we find from (30) and 
(35) that the squared coherence for LSFR and RSFR decreases inverse proportionally to N p and N p , 
respectively. This implies that more pulses are required by RSFR to achieve the desired performance 
with all other parameters, i.e., M t , N r and M, being equal. 

Before ending this section, we note that the expectation of the approximation error in (24) can be 
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represented by 

(36) 

where one can see that a decrease in the product of aw,--, (seen in (31)) and the radian frequency step, 
2nAf, increases both the approximation error and the squared coherence. Given aS,. v , an increase in 
A / would reduce the approximation error pw . However, this would increase the bandwidth required by 
RSFR. 



IV. Decoupled estimation of angle, velocity and range with reduced complexity 

Solving the l\ minimization problem of (11) requires polynomial time in the dimension of s. For the 
discretization discussed in Section II, the joint estimation of angle, velocity and range requires complexity 
of 0((N a NbN c ) 3 ) [3 1] [28]. For large values of N a , Nb and N c , the computational cost of the CS approach 
would be prohibitive. In the following, we propose a decoupled angle-velocity-range estimation approach 
which reduces the search space and thus the computational complexity. 

The scheme needs some initial rough estimates of angle and range. One way to obtain those estimates 
is to use the MFM, which requires forwarding to the fusion center Nyquist sampled data from one pulse. 
In the following, all N r nodes in the system sample all received pulses in a compressive fashion, except N r 
nodes, which sample the first received pulse at the Nyquist rate and all remaining pulses in a compressive 
fashion. Those Nyquist rate samples will be used to obtain coarse estimates of angle and range via the 
MFM. 

The fusion center performs the following operations (also see Fig. 1). 

1) STEP1: Angle and range estimation 

This step uses the first pulse forwarded by each receive node. A fine grid, (a„, , c„, ), . . . , (a„ Ki , c„ x ), 
is constructed around the MFM initial estimates. Then the sensing matrix is constructed as 



0! = Q„ [e^M-** HK, »/c Ci ^ Xv^, ),..., e j2nM - 2c ^ + "< (fl "*. ))/c C Xv m (a njr , )] (37) 
where 



I "1 I 1 | A i . 

L cT s J L cTr J 



I L+L , I — 1, . . . ,N r , m = \ 
<&lm = { „ ■ (38) 

the measurement matrix of size M X (L + L), otherwise 



January 17, 2011 



DRAFT 



14 



The received signals, ?n, . . . ,fjv r i, are stacked in a vector, i.e., 

fi = 0is + ni (39) 

where f i = [ff,, . . . ,f]T l ] T . By applying the Dantzig selector to (39), new and refined angle-range 
information is obtained. 

Thanks to the initial estimates, the search area in the angle -range plane is significantly reduced and 
thus the computational load of CS is lightened. Due to the fact that only one pulse from each receive 
node is used, the range resolution at this step is limited by ^, where B is the signal bandwidth. 
The obtained range estimates will be refined in the next step in which the fusion center will jointly 
process the entire pulse train. Also, due to assumption (Al), Doppler information cannot be extracted 
at this step. 

2) STEP 2: Range resolution improvement and Doppler estimation 

In this step the fusion center processes the entire pulse train forwarded by each receive node. The 
range space around the range estimates obtained in Step 1 is discretized into finer grid points. 
Based on a discretization of the Doppler space, the refined range grid points and the angle estimates 
obtained in Step 1, i.e., (a mi ,b mi ,c mi ), . . . ,(a mK ^b mKi ,c mK ^), the fusion center formulates a sensing 
matrix and extracts angle-Doppler-range information in a CS fashion. 

To further reduce the complexity of CS reconstruction, the MFM can be applied before CS to provide 
angle-Doppler-range estimates around which a finer grid can be constructed and used by CS. In that 
case MFM would be applied based on the grid points (a m , b m , c mi ), . . . , (a mK , b mK , c mK ). 
For the case in which there are stationary targets and moving targets, the angle estimation can be 
further improved by taking into account Doppler information. 
Assuming that the MFM is used for initial estimation, the complexity of two steps is respectively 
0(N a N c (N r L+(N r -N r )M)+Kl) and 0(K 2 (N r (L-M)+N r N p M)+Kl), where K 3 is the number grid points 
used by CS at Step 2. Generally, it holds that K\+K] « N a N c {N r L+(N r -N,)M)+K 2 (N r (L-M)+N r N p M) 
for a small number of targets. Therefore, the computational load is mostly due to the initial estimation. As 
compared to the complexity of the joint angle-Doppler-range CS approach, i.e., 0((N a NbN c ) 3 , considerable 
computations can be saved by using the proposed decoupled scheme for large values of N a ,Nb and N c . 

The computation savings, however, may be obtained at the expense of detection accuracy, unless the 
initial estimates provided by the initial estimation method are reliable. Reliable estimates here refer to 
the initial estimates whose distances to the true target locations are within the resolution cell that is 
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determined by the initial estimation. Then all the targets can be captured based on the finer angle -range 
grid points constructed around the reliable initial estimates. For the instance of the MFM, the performance 
in providing good initial estimates depends on several factors; (i) the signal-to-interference ratio (SIR), 
which can be improved by employing more data; (ii) angular, range or Doppler resolution, which is 
improved by increasing N r or N p ; (iii) the distance between the adjacent grid points. (In the worst case in 
which the targets fall midway between grid points, the targets may fail to be captured by the closest grid 
points if the spacing of adjacent grid points is too large. An empirical approach to select grid spacing 
was discussed in [19]. That approach is also applicable to the MFM); and (iv) the threshold for hard 
detection. A small threshold should be used in order to reduce the miss probability. However, this implies 
that more grid points need to be considered for the CS approach following the MFM as compared to 
a larger threshold. In summary, the performance of the MFM can be improved at the expense of more 
transmit power and increased complexity. 

V. Simulation Results 

We consider a MIMO radar system with transmit and receive nodes uniformly distributed on a disk of 
radius 10m. The carrier frequency is / = 5GHz. Each transmit node uses orthogonal Hadamard waveforms 
of length L = 512 and unit power. The received signal is corrupted by zero-mean Gaussian noise. The 
signal-to-noise ratio (SNR) is defined as the inverse of the power of thermal noise at a receive node. A 
jammer is located at angle 7° and transmits an unknown Gaussian random waveform. The targets are 
assumed to fall on the grid points. Throughout this section, the CS approach uses a measurement matrix 
with Gaussian entries. 

A. Range resolution of the CS-based SFR and conventional SFR 

In this subsection we provide some simulation results to show the superiority of CSSF MIMO radar as 
compared to MFSF MIMO radar in terms of range resolution. Figure 2 shows the normalized amplitude 
estimates of target reflection coefficients for CSSF MIMO radar and MFSF MIMO radar in one realization. 
Since the multiple colocated antennas fail to improve range resolution, we consider a single transmit and 
receive antenna here for simplicity. Let M = 10, N p = 30 and the carrier frequencies be randomly selected 
within the frequency band [5, 5.029]GHz. The CSSF radar uses 10 measurements per pulse while MFSF 
radar obtains 665 measurements per pulse. Various values of SNR are considered. The spacing between 
two adjacent grid points is 2m. There are six targets at ranges [1024, 1028, 1032, 1036, 1040, 1044]m. 
Figure 2 shows that the peaks corresponding to all targets can be distinguished from each other for the 



January 17, 2011 



DRAFT 



16 



CSSF radar while for the MFSF radar some peaks are lost. This verifies the observations of Section III-A 
that CSSF radar has the potential to achieve higher range resolution than does MFSF radar. 

B. Range estimation for CSSF MIMO radar 

The goal of this subsection is to test the performance of CSSF MIMO radar based on LSF and RSF 
Figure 3 compares the numerical and theoretical squared coherence of the sensing matrix corresponding 
to two adjacent grid points in the range plane for different numbers of pulses and various values of the 
linear frequency step Af = I MHz, 4MHz and SMHz. All the results shown in Fig. 3 are the numerical 
squared coherence averaged over 100 independent and random runs and the theoretical squared coherence 
for LSFR and RSFR calculated based on (30) and (34). We consider the case in which M, = M = 10, 
N r = 1 and the grid step is Ac = 7.5m. For a fair comparison, we choose random step frequencies within 
the same frequency band as in LSFR, i.e., / + [0, (N p - 1)A/]. It can be easily seen that the numerical 
squared coherence of the sensing matrix for LSFR perfectly matches with the theoretical results in (30). 
The numerical squared coherence of the sensing matrix for RSFR approaches the theoretical results in 
(34) as Af increases and approaches 1 /N p as the number of pulses increases. It is also verified by Fig. 
3 that LSFR exhibits lower coherence of the sensing matrix than does RSFR. 

Figure 4 shows the receiver operating characteristic (ROC) curves of the range estimates produced 
by the random and linear step-frequency technique based on 200 random and independent runs. Here, 
the probability of detection (PD) is the percentage of cases in which all the targets are detected. The 
probability of false alarm (PFA) is the percentage of cases in which false targets are detected. We consider 
a case in which the angle and speed of three targets are the same and assumed to be known. In each 
independent run, the target angle and speed are randomly generated. The ranges of three targets are fixed 
to 1005m, 1010m and 1045m, respectively. The power of the jammer signal is 4 and SNR= OdB. We can 
see that the use of LSF yields better performance than randomly choosing the carrier frequency within 
the same frequency band. In this particular case, CS-based RSFR requires 12 pulses to generate the ROC 
performance that can be achieved by CS-based LSFR using only 9 pulses. The performance of LSFR 
and RSFR based on the MFM is also shown in Fig. 4. It can be seen that the former using 12 pulses 
is far better than the latter with the same number of pulses. It can also be seen that CSSF MIMO radar 
outperforms MFSF MIMO radar. 
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C. The joint angle-Doppler-range estimation of CSSF MIMO radar 

Figure 5 shows the ROC curves of the angle-speed-range estimates yielded by CSSF MIMO radar 
using the decoupled scheme. The angle-speed-range estimates have been obtained based on 200 random 
and independent runs. The cases in which M, = 10, N r - N r - 7 and N p = 12 are shown in Fig. 5. 
The azimuth angle and range of three targets are randomly generated in each run but the spacing of 
angle and range between targets are fixed to 0.3° and 7.5m, respectively. The speeds of three targets are 
lOm/s, 30m/ s, and 60m/ s. The power of the jammer signal is 4 and SNR= OdB. The performance of 
MFSF MIMO radar, shown in Fig. 5, is obtained in the same decoupled fashion, i.e., 1) estimate target 
angle and range based on a single pulse; then refine the angle estimates based on the finer angle grid 
points around the initial angle estimates by using the MFM; and then 2) process the entire pulse train 
to extract angle-speed-range information by discretizing the speed space, constructing finer range grid 
points around the initial range estimates and utilizing the initial angle estimates obtained in 1). One can 
see that MFSF MIMO radar is inferior to CSSF MIMO radar even when using far more measurements 
than the latter. 

VI. Conclusions 

We have presented a CSSF MIMO radar system that applies SF to CS-based MIMO radar. The technique 
of SF can significantly improve range resolution. We have shown that CSSF MIMO radar has the potential 
to achieve better resolution than MFSF MIMO radar, and that more pulses are required by RSFR than by 
LSFR to achieve the desired performance with all other parameters being the same. The angle-Doppler- 
range estimation requires discretization of the angle-Doppler-range space into a large number of grid 
points, which would increase the complexity of the CS approach. We have presented a CSSF MIMO 
radar scheme that by decoupling angle -range estimation and Doppler estimation achieves significant 
complexity reduction. The proposed technique applies to slowly moving targets and relies on initial 
rough angle-range estimates. Assuming that the initial estimates do not miss any targets, the proposed 
low complexity scheme maintains the high resolution of the CS approach. 
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Fig. 1. Schematic diagram of the proposed decoupled scheme. 




Fig. 2. Normalized amplitude estimates of target reflection coefficients for the CSSF radar and MFSF radar (one realization 
for M, = N r =l and N p = 30). 
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Fig. 3. Average squared coherence of the sensing matrix for different numbers of pulses N p over 100 independent and random 
runs (A/ = 4 x 10 6 , M = M, = 10 and N r = 1). The distance of two grid points in the range plane is Ac = 7.5m. 




Fig. 4. ROC of range estimates obtained with linearly and randomly stepped frequency CSSF MIMO radar and MFSF MIMO 
radar (M = M, = 10, N r = 1 and A/ = 1MHz). 
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Fig. 5. ROC of target detection based on angle-speed-range estimates yielded by the proposed decoupled scheme in Section 
IV for CSSF MIMO radar and MFSF MIMO radar (M, = 10, N r =N r = l and N p = 12). 
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